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ABSTRACT 

During the September-October 2008 outburst of the accreting millisecond pulsar SAX J1808. 4-3658, the 
source was observed by both Suzaku and XMM-Newton approximately 1 day apart. Spectral analysis reveals a 
broad relativistic Fe Ka emission line which is present in both data-sets, as has recently been reported for other 
neutron star low-mass X-ray binaries. The properties of the Fe K line observed during each observation are 
very similar. From modeling the Fe line, we determine the inner accretion disk radius to be 13.2 ±2.5 GM/c 2 . 
The inner disk radius measured from the Fe K line suggests that the accretion disk is not very receded in the 
hard state. If the inner disk (as measured by the Fe line) is truncated at the magnetospheric radius this implies 
a magnetic field strength of ~ 3 x 10 s G at the magnetic poles, consistent with other independent estimates. 

Subject headings: stars: neutron — X-rays: binaries — X-rays: individual: SAX J1808. 4-3658 
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1. INTRODUCTION 

The accreting millisecond pulsar SAX J1808. 4-3658 was 
the first X-ray binary discovered to pulsate i n the mil- 
lisecond range dWijnands & van der Klisl 1 1998). It is a 
transient compact binary wit h a ~2.01 hr orbital period 
(Chakrabartv & Morgan 1998) and an outburst recurrence 
time of approximately 2.5 years. SAX J1808. 4-3658 has 
been seen in outburst six times since 1996, the most recent 
being in September-October 2008. The presence of pulsa- 
tions implies the accreting gas is channeled onto the magnetic 
poles producing hot spots and accretion shocks. The rotation 
of the neutron star modulates this emission producing pulsa- 
tions that add to the impuls ed accretion disk and Compt onized 
emission (see for example [Poutanen & Gier lihski 2003|) . 

A tool to examine SAX J1808.4-3658 is the presence of 
an iron K-shell fluorescence emission line in its X-ray spec- 
trum, similar to th ose observed in many other X- ray bina- 
ries and AGNs (see lMiUdl2007t iFabian et alJl2000l for com- 
prehensive reviews). It is generally thought that the shape 
of iron-K emission line formed in the inner region of ac- 
cretion disks around black holes (both stellar-mass and in 
AGN) and neutron stars is sensitive to the distance of the 
line forming region from the compact object - as this re- 
gion gets closer to the compact object relativistic Doppler 
effects and gravitational re dshifts become str onger, produc- 
ing an asymmetric profile (Fabian et al] |1989l) . The iron K 
emission line can therefore provide an important way to ex- 
amine the inner accretion disk. While this line has lon g been 
used to study the inner disk around black holes (e.g., Miller 
120071) . only recently has it been shown that the lines in neu- 
tron star low-mass X-ray binaries also display the characteris- 
tic asymmetric shape too ( Bhattach aryya & S trohmay erl2007l: 
ICackett eT al. 2008b). The detection of an iron line in an ac- 
creting millisecond pulsar is useful for constraining the mag- 
netospheric radius and hence the strength of the magnetic field 
in this neutron star. 

During the Sept/Oct 2008 outburst, SAX J1808.4-3658 
was observed by both Suzaku and XMM-Newton, approxi- 
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mately 1 day apart. Initial brief rep orts of the F e K lin e 
seen by both telescop es were given by Cac kett et al.l (|2008a); 
IPapitto et ail d2008al) . Furt her analysis of the XM M-Newton 
data has been presented by IPapitto et al.l d2008bl) . Here, we 
present a joint analysis of both Suzaku and XMM-Newton data 
that yields the tightest constraints on the Fe K emission line 
present in the spectra. 

2. DATA REDUCTION AND ANALYSIS 

The light curve of the Sept/Oct 2008 outburst is given in 
Fig. [U The daily average count rates are taken from the 
SwiftfBAT Hard X-ray Transient Monitor in the 15-50 keV 
energy range. The times when the Suzaku and XMM-Newton 
observations occurred are also indicated. Below we describe 
the data reduction for both the Suzaku and XMM-Newton, be- 
fore detailing our spectral analysis. 

2.1. Suzaku data 

Suzaku observed SAX J1808.4-3658 starting at 2008 
October 2, 16:32 UT, for approximately 42.5 ks (ObsID 
903003010). The XIS was operated in 1/4 window mode, 
with a 1-sec burst option, and the telescope was pointed at 
the nominal XIS position. The XIS detectors, which cover 
the energy range from approximately 0.5-10 keV, were oper- 
ated in both 3x3 and 5x5 edit modes during the observation. 
For each working XIS unit (0, 1 & 3), we extracted spectra 
from the cleaned event files using xselect. For the front- 
illuminated detectors (XIS and 3), the extraction region used 
was a box of size 270" by 330" centered on the source. We 
choose to use a box rather than a circular extraction region 
as the used region of the detector is a thin 256 by 1024 pixel 
strip (approximately 270" by 1080") when using 1/4 window 
mode. The background was extracted from a 250" by 230" 
box situated towards the end of the 1/4 window so as to be 
free from source photons. The response files were generated 
using the xisresp script which uses the xisrmfgen and 
xissimarfgen tools, xissimarfgen is used with no 
binning of the response file and simulating 200000 photons. 
We then added the spectra (and averaged the responses) from 
the two front illuminated detectors (XIS and 3) together us- 
ing the addascaspec tool, as recommended by the Suzaku 
XIS team. The summed good exposure time of the XIS + 
3 spectrum is 42.4 ks (this takes into account the 1-sec burst 
option which leads to a 50% livetime fraction). 
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FIG. 1.— Swift/BAT 15-50 keV light curve of the September-October 
2008 outburst of SAX J1808.4-3658. The dashed and solid lines mark the 
mid-points of the XMM-Newton and Suzaku observations respectively. 

In this Letter, we choose not use data from the back- 
illuminated XIS 1 detector. As the back-illumination in- 
creases the effective area at soft energies, this detector is more 
prone to pile-up. We found that pile-up was present in the 
XIS 1 observation (but not significantly in the other detec- 
tors). We tested using an annulus source extraction region that 
avoids the piled-up central region of the source, but found that 
when we did this the signal-to-noise ratio was significantly 
decreased, especially through the Fe K band (6.4-6.97 keV) 
of interest here. 

In addition to the XIS spectra, we also extract the spec- 
trum from the PIN camera which is part of the Hard X-ray 
Detector (HXD). We extracted the PIN spectrum from the 
cleaned event file following the standard analysis threads on 
the Suzaku website. The PIN non-X-ray background was ex- 
tracted from the observation-specific model provided by the 
instrument team, and was combined with the standard model 
for the cosmic X-ray background. 

2.2. XMM-Newton data 

XMM-Newton observed SAX J1808.4-3658 starting on 
2008 September 30, 23:53 UT, for a total of -63 ks (Ob- 
sID 0560180601). Here we present data from the PN detec- 
tor, which was operated in timing mode during the observa- 
tion to prevent pile-up. Processing of the ODF failed using 
epproc, therefore we used the PN event file provided by the 
XMM-Newton team. We checked for background flares, and 
found that the first —5 ks were affected by a flare, therefore 
we excluded that section of the data from the rest of the analy- 
sis. We extracted the spectrum using xmmselect. In timing 
mode the detector is continously read-out in the Y direction, 
leading to a bright streak on the detector. Thus, the source ex- 
traction region was chosen to be of width 24 pixels centered 
on the source and with the length covering the entire readout 
streak. The background spectrum was extracted from a strip 
at the edge of the observed window. The good exposure time 
of the source spectrum was 43.4 ks. The response files were 
generated using the arf gen and rmf gen tools. 

2.3. Spectral Fitting 

We fit the spectra using XSPEC version 1 1 (lArnaudll 19961) . 
First, we fit both the Suzaku and XMM-Newton spectra sep- 
arately. For the Suzaku data we fit the XIS spectrum over 
the 0.7 - 10 keV energy range, and the PIN spectrum over 
the 14-45 keV range. For the continuum model, we used 
an absorbed multicolor disk-blackbody + a single tempera- 
ture blackbody and an additional power-law (e.g. iLin et al.1 



Parameter 


Suzaku 


XMM 


N H (10 21 cm 2 ) 


0.46 ±0.06 


2.3 ±0.1 


Disk T in (keV) 


0.48 ±0.01 


0.23 ±0.01 


Disk normalization 


590 ±51 


20400^™° 


Blackbody, kT (keV) 


1.00 ±0.04 


0.40 ±0.01 


Blackbody normalization (10 -3 ) 


1.27 ±0.02 


1.25 ±0.06 


Power-law index 


1.93 ±0.02 


2.08 ±0.01 


Power-law normalization 


0.20 ±0.01 


0.33 ±0.01 


Enne (keV) 


640+ o.06 


640+ o.06 


13 


-2.9 ±0.3 


-3.0 ±0.2 


in 


50i" 
12.75' 3 


59±4 


Rin (GM/C 2 ) 


13.0=1=3.8 


EW (eV) 


134 ±30 


1 18 ± 10 




1.17 (2550) 


1.09 (1950) 



NOTE. — All uncertainties are at the 90% confidence level. The Suzaku 
spectra are fit over the range 0.7-45 keV, and the XMM-Newton spectra are fit 
from 1.2-11 keV. 



TABLE 2 

Iron line parameters from joint fitting 



Parameter Suzaku & XMM 



E Une (keV) 6.40+ 003 
(3 -3.05 ±0.21 

i (°) 55% 
S, n (GM/c 2 ) 13.2 ±2.5 

EW, Suzaku (eV) 144+^ 
EW,XMM(eV) 113 ±10 

X]M) 1-14(4502) 

NOTE. — All uncertainties are at the 90% confidence level 

l2007t ICackett et al.ll2008bh . All components are required by 
the data and statistically improve the fit. For the absorption 
we used the phabs model. In the Suzaku/XIS spectrum we 
found that there is a detector feature at approximately 1.85 
keV. We attribute this to a change in the instrument response 
that is not yet corrected for by the response files (similar fea- 
tures havejje^nse^njn other 1/4 window mo de observations, 
e.g.. iMiller eTail 120081: ICackett et al.l l2008b). We therefore 
model the feature with a Gaussian absorption line so that it 
does not effect the continuum fit. We allowed a constant be- 
tween the Suzaku/XIS and Suzaku/PIN spectra to account for 
any offset in flux calibration, finding c= 1.37 ±0.05. Such 
an offset between XIS and PIN flux can occur in 1/4 window 
mode because small extraction regions must be used. This is 
a known problem detailed in the 'Calibration Issues' on the 
Suzaku website. 

Initially, we fit the continuum ignoring the iron line region 
4-7 keV. When examining this region after having fit the 
continuum, an i ron line was clearly present, thus we added 
the diskline (Fa bian et al.| [l989) model to fit the line and 
re-fit all the parameters. The energy of the line center was 
constrained to be within the Fe K band range, 6.4 to 6.97 keV. 
We left the power-law disk emissivity profile, j3, and inner 
accretion disk radius, Rj„, as free parameters, while the outer 
accretion disk radius was fixed at a large value (1000 GM/c 2 ). 
The inclination, i, of the disk was also a free parameter, yet 
was constrained to lie within 36 to 67 degrees - the possi- 
ble range of inclination s determined by op tical observations 
of SAX J1808.4-3658 (iDeloye et alJl2008h . The best-fitting 
parameters (continuum & line) are given in Tab. Q] 

We used the same method when fitting the XMM- 




TABLE 1 

Spectral parameters from separate fitting 



Broad Iron Line in SAX J1808.4-3658 
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FIG. 2.— Suzaku/XlS 1+3 (black), Suzaku/PJN (red) and XMM-Newton/PN 
(green) spectra of SAX J1808. 4-3658 (top panel). The residuals of the best 
fit (x = (data-model)/<r) are plotted in the bottom panel. 
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FIG. 3. — Iron K emission line in SAX J1808.4-3658 detected by Suzaku 
(black) and XMM-Newton (red). The ratio of data to continuum model is 
plotted to show the iron line profile. The black, solid line shows the best- 
fitting iron line model to the Suzaku data, and the red, dashed line shows the 
line model for the XMM-Newton data. 



Newton/PN spectrum. When examining the PN spectrum we 
found significant residuals below 1.2 keV, therefore, we only 
fit the spectrum over the 1.2-11 ke V energy range. As with the 
Suzaku data, a clear Fe K line was present in the data, which 
we again fitted using the diskline model. The best-fitting 
parameters (continuum & line) are given in Tab.Q] and the Fe 
K line profiles are shown in Fig. [3] 

The continuum shapes from Suzaku and XMM-Newton are 
different. This is not unexpected given that the observations 
were non-simultaneous, and the flux (see Fig. Q]) had changed 
between the observations. However, the Fe K line properties 
are similar (all the parameters are consistent within the uncer- 
tainties), thus, we chose to fit the Suzaku and XMM-Newton 
spectra simultaneously to place tighter constraints on the inner 
disk radius. In this joint fit, we allowed the continuum param- 
eters to be completely independent, yet tied all the line pa- 
rameters (except the normalization) between the two datasets. 
Parameters from this joint fit are given in Tab. [2] We find an 
inner accretion disk radius of Rj„ = 13.2 ±2.5 GM/c 2 . The 
inclination of the disk is measured to be i = 55+ 8 degrees. We 
show the broadband spectrum in Fig. [2] 

In addition to this phenomenological approach, we also in- 
vestigated a number of self-consistent spectral fits, including 
a relativistic disk reflection spectrum. The Suzaku spectra can 
be fit well using the ref sch model instead of a power-law 
component, leading to a reflection fraction of 0.5 ± 0. 1 . How- 
ever, we note that there is little evidence of the Compton back- 



scattering hump expected from disk reflection; the HXD/PIN 
spectrum can be fit well using only a simple power-law. We 
will investigate disk reflection in neutron star systems in an 
upcoming paper. 

2.4. Estimating the magnetic field strength 

If we assume that the inner accretion disk is truncated at 
the magnetospheric (Alfven) radius (where magnetic pressure 
balances the ram pressure from infalling material), then from 
the source luminosity and reasonable assumptions about mass 
and radius we can esti mate the magnetic field strength (for 
example see Eq 6 .19 in | Frank et al.ll2002h. Her e we modify 
the formulation of Ibragimov & Poutanen (2009) substituting 
Ri n = xGM/c 2 into their Eq. 14 to give the following expres- 
sion for the magnetic dipole moment: 



fi = 3.5xl0 23 k- 7/4 x y 4 

fang Fbol 



M 



1.4 M 







rj 10 9 erg cm 2 s 



1/2 



D 



3.5 kpc 



Gem 3 (1) 



where r/ is the accretion efficiency in the Schwarzschild 
metric, f ang is the anisotropy correction factor ( which is 
close to unity, see Ibragi mov & Poutanen! 120091 for de- 
tails), and the coefficient &a depends on the conversion 
from spherical to dis k accretion . Nu merical simulations 
suggesting fc^ = 0.5 dLpng et al.l l2005h . and as noted by 
Ibragimov & Poutanen (2009) , theoretical models predict 
k A < 1.1 dPsaltis & Chakrabartvll999HKluzniak & Rappaporil 
2007). For ease of comparison with previous magnetic 
field estimates, we assume a distance D = 3.5 ±0.1 kpc 
dGallowav & Cumming1l2006|) and mass M = 1.4 M Q . 

We also use a bolometric conversion fa ctor of 
Fbol/Fi-zSkeV = 2.12 dGallowav & Cummin d 120061) . The 
2-25 keV fluxes that we measure from the Suzaku and 
XMM-Newton observations are (1.12 ± 0.01) x 10" 9 erg 
cm" 2 s" 1 and (1 . 16 ± 0.01) x 10~ 9 erg cm" 2 s" 1 , respectively. 
Here, we just use the average of the two, which leads to 
F bo i = (2.42 ± 0.02) x 10" 9 erg cm" 2 s" 1 . 

Using Rj„ from the joint SuzakulXMM-Newton Fe line fit, 
along with assuming k& = 1, fang = 1> an d rj = 0.1 leads to 
/U = (1.6 ±0.5) x 10 26 G cm 3 , taking into account the uncer- 
tainty in the inner disk radius, flux, and distance. For a stellar 
radius, R = 10 km, this leads to a magnetic field strength of 
B = (3.2 ± 1 .0) x 10 8 G at the magnetic poles. Note if we as- 
sume k A = 0.5, we would infer a larger magnetic field strength 
approximately a factor of 3 larger. 

3. DISCUSSION 

We have observed a broad, relativistic Fe K emission line 
in the accreting millisecond X-ray pulsar SAX J1808.4-3658, 
using both Suzaku and XMM-Newton. The line is similar in 
both observations, and by modelling the iron line we deter- 
mine the inner disk radius, R,„ = 13.2 ±2.5 GM/c 2 . Mak- 
ing reasonable assumptions about the neutron star radius and 
mass allowed us to estimate the magnetic field strength B = 
(3.2 ± 1.0) x 10 8 G at the magnetic poles. The broadband 
spectrum from Suzaku displays a hard power-law tail that ex- 
tends to at least 45 keV. Although we were able to fit the 
broadband spectrum well with a reflection model, we note 
there is no clear sign of the Compton hump that can be seen 
in black hole systems. Neutron star LMXB spectra are cer- 
tainly not as simple as black hole LMXB spectra, with the 
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former often requiring multiple thermal components, which 
may lead to hide the Compton hump. 

In comparison with our inner disk rad i us, from fitting the 
XMM-Newton data alone, iPapitto et all d2()08h) find R in = 
8.7^27 GM/c 2 consistent with our value, though a little 
smaller. The small difference may be attributed to the fact 
that while we fixed the outer disk radiu s at a large value (so 
that is not does affect the line profile), Papi tto et al.l {2008b) 
left it as a free parameter. If we also leave the outer disk ra- 
dius as a free parameter, we reproduce their result. Of all the 
parameters of the line model, the line profile is least sensitive 
to the outer disk radius, R„„ t , and this parameter cannot be 
well constrained from spectral fitting - for a typical emissiv- 
ity profile of 7?" 3 , there is a strong weighting to the innermost 
part of the disk. Thus, the slightly different assumptions about 
the outer disk have only a small effect on the measured inner 
disk radius. 

An indep endent estimate of the i n ner a ccretion disk radius 
is made by llbragimov & Poutanen (2009). From RXTE ob- 
servations of the 2002 outburst of SAX J1808.4-3658, these 
authors study the evolution of the pulse profiles, finding a sec- 
ondary maximum during the late stages of the outburst. They 
attribute this to the accretion disk receding and allowing a 
view of the lower magnetic pole. Under this assumption, they 
estimate the inner disk radius, at that specific time in the out- 
burst, to be 19.5 km for M = 1 .4 M , R = 12 km, and i = 50°, 
with larger estimates for larger inclinations. This is close to 
o ur determinat i on of t he inner disk radius. 

Pand el et"ai1 (|2008;) find a broad, relativistic Fe K line in 
the spectrum of 4U 1636-536. They suggest that the line is a 
blend of several Ka lines in different ionization states. How- 
ever, from their modeling of the line profile with 2 ionization 
state lines (one at 6.4 keV, and one at 7 keV) the inner disk 
radii they determine for each ionization state are very similar 
(within 20 GM/c 2 of each other). Yet, it would seem difficult 
to produce lines with such different ionization states so close 
to each other. In addition, the disk emissivity heavily weights 
the emission from the innermost part of the disk, where the 
highest ionization state line would come from. Thus, even if 
multiple lines are present, it seems reasonable to assume that 
the highest ionization state line would completely dominate. 
Here, we find we do not require a second ionization state to 
model the line profile. 

The power-law index measured during the Suzaku and 
XMM-Newto n observation s, as we ll as the power spectrum 
presented in Papitt oet alj ((2008b), is typical of the low- 
hard state (often referred to as the 'island' state, e.g., see 
Ivan der Klisll2006l for further details on neutron star LMXB 
states). The inner disk radius that we measure from the Fe 
line is only ^2 times the innermost stable circular orbit (for a 
Schwarzschild metric), suggesting that there is little recession 
of the accretion disk in this hard state. This is in di sagree- 
ment with the model proposed by iDone et ail (|2007) where 
the disk is recessed in the island state. A similar result has 
been seen in black hole LMXBs - observations of iron lines 
and/or a cool disk in the low-hard state of black hole LMXBs 
also suggests th at there is little rec ession of the ac cretion disk 
in such a state (Miller et al. 2006a b). Moreover, Rykoff et al. 



d2007l) show that there is a disk component in the black hole 
LMXB XTE J1817-330 that cools as L x ~ T 4 as the source 
transitions to the low-hard state, with an inner disk radius con- 
sistent with the innermost stable orbit at all times. 

Our estimate of the magnetic field strength is broadly con- 
sistent with previous estim ates. From the long-term spin 
down. lHartman et alJd2008l) determine B < 1 .5 x 10 8 G at the 
magnetic poles. iDi Salvo & Burderil d2003l) combine the qui- 
escent luminosity of SAX J1808. 4-3658 along with consider- 
ations of the magnetospheric radius during quiescence to esti- 
mate B = 1 - 5 x 10 8 G (at the equator, note that for a dipole the 
field strength at the equator is a factor of 2 smaller tha n at the 
poles). An earlier estimate by Psaltis & Chakra bartyl d 19991) 
is less constraining implying a large possible range in surface 
dipole field of 10 8 - 10 9 G at the stellar equator. Moreover, 
llbragimov & Poutanen! (120091) estimate B = (0.4-1.2) x 10 8 
G from their measurement of the inner disk radius. The esti- 
mation of the magnetic field strength here is, of course, depen- 
dent on a number of assumptions. In particular, the coefficient 
Ica is the most uncertain, with a lower value of Ua increasing 
our estimate of the magnetic field strength. The uncertainty 
from our particular choice of k,A is not included in the quoted 
confidence limit. Moreover, our estimate is based on the as- 
sumption that the disk is truncated at the magnetospheric ra- 
dius. This, however, may not necessarily be the case, as MHD 
simulations have shown that periodicities at the star's rota- 
tion frequency can still be pr esent during unstable accretion 
dKulkarni & Rom anova 2009). Additionally, the presence of 
the stellar magnetic field affects the inner part of the accretion 
disk which may in turn affect the iron line profile. It is not 
clear how this will change the shape of the line profile from 
the standard diskline model used here. 

Finally, optical observations during the 2008 outburst of 
SAX J1808. 4-3658 have lead to an improved mass function 
for th e pulsar, with M\ sin 3 / = 0.48!q;}4 M q dElebert et alJ 
2009). From these observations they are unable to constrain 
the binary inclination any better than previously. If we as- 
sume that the inclination we measure from the Fe K line is the 
same as the inclination of the binary system, then this would 
imply quite a low pulsar mass of 0.9^ 4 M©. However, there 
is still a large uncertainty in this mass determination, which is 
consistent with the canonical 1.4 Mq. In contrast, tight con- 
straints on the thermal emission from the neutron star surface 
during quiescence imply that SAX J1808.4-3658 requires en- 
hanced levels of core coolin g, which can poss ibly be achieved 
through have a larger mass dHeinke et al . 2008). 

It is also interesting to note that Eleb ert et alj Q2009) esti- 
mate a shorter distance to SAX J1808.4-3658 of ~ 2.5 kpc 
based on the equivalent width of the interstellar absorption 
lines. This shorter dis tance is the same as tha t estimated from 
type-IX-ray bursts by lin 't Zand et a l. (2001). Such a distance 
would reduce our magnetic field strength estimate by 0.7. 
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